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Abstract. We prove a multidimensional extension of Selberg's central 
limit theorem for the logarithm of the Riemann zeta function on the 
critical line. The limit is a totally disordered process, whose coordinates 
are all independent and Gaussian. 



1. Introduction 

A classical result of Selberg [10] (see also Laurincikas, [6]) states that 
the classical continuous determination of the logarithm of the Riemann zeta 
function is asymptotically normally distributed, in the sense that if T is a 
regular Borel measurable subset of C, 

-2T 



TJ t ./^ log log X J 2 *Jr 



lim 

T->oo T 

where l is the indicator function, and regular means that the boundary of 
r has zero Lebesgue measure. 

If we let 

logC^ + iue"*) 

L\(N,u) := -== 

VlogiV 

then Selberg's result implies that 

lim / i{L A (iv,tt)Gr} du = p{G A er} 

where G\ = g{ 1} +iGf ] is a complex- valued Gaussian random variable with 
mean zero and variance A/2, i.e.: G x and G{ are independent, centered, 
andE[(Gi 1) ) 2 ]=E[(Gi 2) ) 2 ] = A/2. " 

It is now a natural question, at least from a probabilistic standpoint, to 
look for an asymptotic distribution for (L\ 1 (N, •),..., L\ k (N, •)), for differ- 
ent A,;'S. 
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> Xk > 0, and for every (Ti,i < k) regular, 

k 

, L Xk (N, u) g r fc } dn = n P { G\ J eTj}. 

i=i 

(1) 

We now note that if (p x = d[ 1] + \d[ 2 \ A > o) is a totally disordered 
complex- valued Gaussian process, meaning that ^D^\ A > 0^ and 

, A > 0^ are two independent Gaussian processes all of whose coordi- 
nates are independent with E[(D^ 1) ) 2 ] = E[(d[ 2) ) 2 ] = A/2, then the quantity 
on the right hand side of (P) is 

¥{D Xl 6 Ti, . . . , D\ k £T k }. 

Theorem Q] is an attempt to move from the deterministic set up of the 
Riemann zeta function, and the "static" central limit theorem of Selberg 
into a more "dynamic" probabilistic world, where a process appears in the 
limit. However, this process is quite wild. In the next section, we comment 
about it, and some of its occurrences in random matrix theory. Finally, in 
the third section we prove Theorem [T] using the method of moments. 

Remark. Our methods apply equally well to any L-function from the Selberg 
class, but for concreteness and for the sake of simplicity we only state here 
the result for the Riemann zeta function. 

2. Some remarks on total disorder process 

2.1. Non-measurability of the total disorder process. The total dis- 
order process is a "wild" process; indeed there is no measurable process 
(A, a;) i— * D\(lo) which would be a modification of (D\,\ > 0), i.e. ^{D\ = 
D\) = 1 for all A. Indeed, if so, we would get (use Fubini) 

b _ 

I>\dA = a.s., 
hence 

D x = dAdP, 

which is absurd (for some further discussion on the total disorder process, 
see page 37 of [7]). 

2.2. The total disorder process in random matrix theory. The total 
disorder process has already been observed asymptotically in random matrix 
theory, although in a different guise. Let Z v (0) = det(I - Ue-' ld ) be the 



Theorem 1. For X\ > A2 > • • • 

lim / l{L Xl (N,u) GTi,... 
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characteristic polynomial oi an N x N unitary matrix U chosen with Haar 
measure, then Hughes, Keating and O'Connell [4] prove that 

lo g Zu(e) 

i log iV 

weakly converges to X(9) + \Y{9), where X(9),Y(9) are independent Gauss- 
ian processes with covariance structure 



E[X(9 1 )X(9 2 ))=E[Y(9 1 )Y(9 2 )) 



1 if 0i = 6» 2 
otherwise 



This was used to provide an explanation for the covariance structure of 
Cjj(s,t), the number of eigenangles of U that lie in the interval (s,t), found 
earlier by Wieand [121 [13] . A separate explanation was given by Diaconis 
and Evans [2]. Let 

Cu(s,t) - (t-s)N/2ir 



Cu(s,t) :-- 



Wieand proves that for fixed s, t, if the matrices U are chosen with Haar 
measure from the unitary group, then Cu(s,t) converges in distribution, as 
N — > oo, to a standard normal random variable. In fact she goes much fur- 
ther by proving weak convergence of Cu(s,t) to a certain Gaussian process 
C(s,t). 

Theorem 2 (Wieand). For —it < s < t < it, the finite dimensional distri- 
butions of the process Cu(s,t) converge as N — > oo to those of a centered 
Gaussian process C(s,t) with covariance structure 



ifs = s',t = t' 
ifs = t',t = s' 

if s = s' or if t = t' but not both 
if s = t' or if t = s' but not both 
otherwise 

A similar process result had previously been found by Costin and Lebowitz 
[lj for GUE matrices, and Soshnikov [llj considers a process result for count- 
ing the number of eigenangles in an interval with a given minimum displace- 
ment. The surprising thing about these correlations is that they imply that 
even if an interval / contains more than the average number of eigenangles, 
then any subset of / not sharing a common endpoint with I will usually 
still contain its average number. Also, no matter how close two intervals I 
and J are, unless they share an endpoint, then Cj and Cj (with obvious 
notations) are independent. 



E{C{s,t)C(s',t')} 
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Of course, the results of Wieand and of Hughes, Keating and O'Connell 
are strongly related, because 

~ 3mlog Zu(t) _ 3m\ogZu(s) 

3. Proof of theorem Q] 

We will prove Theorem [1] via the method of moments. The following 
lemma will play an essential role in our argument. 



Lemma 3. A complex random variable Z = X + \Y has moments 

\n\2 n a 2n if m = n 
I otherwise 



E [Z m Z r 



(2) 



if and only if X and Y are independent and distributed according to the 
normal law with mean zero and variance a 2 . 



Proof. Let Z = X + \Y where X and Y are independent, centered normal 
random variables with variance a 2 . Consider the joint moment generating 
function of Z and Z: for (a, (3) G C 2 , 



E 



E 



(a+f3)X+i{a-f3)Y 



= e (a+/3) 2 <7 2 /2-( a -/3)V 2 /2 

^— ' n! n! 

n=0 

which is the two- variable moment generating function of ([2]). Conversely, 
assume that ([2]) holds for the joint moments of Z and Z. Then working up 
the above chain of equalities proves that Z = X + iY where X and Y are 
independent, centered gaussians with variance a 2 . □ 

From Lemma El if one can show that for any positive integers k and any 
integers mi, . . . , m^; m, . . . , n^, and if for any Ai > • • • > A& 



E 



-i if 



m t ,n t 



then one may conclude that (D\, A > 0) is a centered complex- valued Gauss- 
ian totally disordered process with covariance structure 

A, if A,- = A 



and E [D Xl D X] ] = for all Ai, Xj. 



otherwise 
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Therefore, Theorem [1] is a consequence of the following 
Theorem 4. Let 

logCq + me^) 

L\(N,u) = j== . 

Vlog N 

If Ai < A2 < • • • < Afc are fixed, then 

/2 k k 
[] L\. (N, u) m iL X3 (N,u) nj dn=H n/.X'/ <){n, r 
- 3=1 3=1 

where S(rrij,nj) = 1 if rrij = rij and zero otherwise. 

We need the following theorem of Selberg, |10| . 

Theorem 5 (Selberg). If n is a positive integer, < a < 1, and T a l n < 
x < T l l n , then 

2n 

1 



rij) 



p<x VP 



dt = (n in e An ) 



r2T 
IT 

for some constant A which depends upon a. 

We also need to calculate the moments of certain prime sums. 
Lemma 6. Given k a positive integer, let Ai > • • • > Afc > 0. Let 

P(X,n) = P(X,n;k,N,u) = , V 



For any non-negative integers mi, . . . , and m, . . . , n^, 

r>2 k k 

m / 

N- 



„2 k k 

lim / TT P(X j ,m j ) m ^P(X j ,n j f j du = TT(n i )! (Xj)^ <5(m,-, nA (3) 

J1 3=1 3=1 



(Note that for the sake of simplicity the variable u does not appear ex- 
plicitly in Q and in some expressions below) 

Proof. We wish to expand out 

\{P(X v m 3 Ti 

3=1 

as a multiple sum over primes. It is exceedingly complicated. We will 
introduce the following notation: For j = 1, . . . , k, let 
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and let 

V j = V (Aj,mj) = jpj : is prime , Pj/ < exp } 

Hence 

k 

[JP(A J -,m j ) m ^ = (logiV)-( mi +- +mfc )/ 2 

exp ( -iu Ylj=i eAfAj logfo.i • • • Pj.mj )) 
X > ^ - 

piePi,...,p fc GP fe yllj=i ll^=i Pj/j 

Similarly, we let 

qj = • • • ><^>j) 

and let 



2j = ^ ( A i> n i) = jqj : is P rime » Qj,e < ex P ( 



40ferij 



and so 

[]P(A i ,n i )" i = (logAf)-( ni+ - +nfc )/ 2 



x 



exp (iu £)J =1 log(9j, i . . . ^ )) 



F( P , q ) = Vex P f^ A lo g( — 



qieSi,...,q fc eS fc \J^j=i 11™/= 1 <7j,£ 

Finally, let p = UjLi Pi an< ^ Q = UjLi Qj ) an d let 

k 

Qj,i ■ ■ ■ Qj,n 3 

c^p l i y -l iug l 

V / U.1 ...I);. 

J'=l 

Therefore, 

[]P(A i ,m i rP(Vi/ J ' = (iogAr)-( m i+-+ m fc+ n i+-+^)/ 2 

3=1 

x £ exp (mF(p,q)) ^ 

, p ,'IS;:::Sia ^/iti- (n^) 

We divide the sum up into two parts, depending on whether F(p, q) equals 
zero or not. The terms where the sum vanishes we call diagonal terms; 
the other terms are off-diagonal The proof of the lemma will follow from 
showing that the off-diagonal terms do not contribute in the large- N limit, 
and using a simple combinatorial enumeration of the diagonal terms, along 
with the prime number theorem, to estimate the diagonal terms. 
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3.1. The diagonal terms. We will see below in §3.21 that since Ai > • • • > 
Afc, then for sufficiently large N, the only way for F(p, q) = is if 



exp ( iV A > lo. 



• • • 1i,n 3 

Pj,l ■ ■ ■ Pj,mj 







for each j = 1 , . . . , k separately. Thus the diagonal terms are those contained 
in the sets 



(5) 



Since pj i and i are both prime, the set T>j is empty unless rrij = rij. Under 
such an assumption, the diagonal terms in (jl]) are 



E 



(pi,qi)eX>i,...,(p fc ,q fc )eX> fc 11? 



IL=i IlZ=i 9j,t 



n e 

3 1 \(Pi.qi)G»; 



1 



■ ■ ■ <ti,nj 



If qj t i, . . . , qj <nj are distinct primes, then there are (rij)\ ways of choosing 
Pj t i, . . . ,Pj t nj such that the products are equal. (If the qj t t are not distinct, 
then the result is similar, but with a different combinatorial factor, and the 
result is at least a couple of logarithms smaller) . Hence 



E 



1 



qj,i ■ ■ ■ qj,n 3 



E 



1 



( 



qj,i ■ ■ ■ qj, nj 



i + o( 



1 



log^iV 



{ nj )\ 



\ 



E 



q<cxp 



^ ! ( l0E (^) +0(I 



(n 



(njV.iXjlogN)^ l + 0( 



l + 0( — j—' 
doe 2 TV' 



l + 0( 



log 2 iV 



1 



dogiv 



(6) 



Hence the diagonal contribution to (Jl|) is 




l + 0( 



1 



loeiV' 



which is the right-hand side of ([3]) in the large- N limit. (The constant 
implicit in the O-term depends on mj,rij,\j and k, but these are all con- 
stants). Hence the proof of the lemma will be complete if we can show 
there is no contribution to ([3]) from the off-diagonal terms, the terms where 
F(p,q) 
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3.2. The off-diagonal terms. Now we show that the non-diagonal terms 
of (|4|) do not contribute to ([3]) in the limit. Upon integrating ([4|) for u 
between 1 and 2, we obtain 

exp (i2F(p, q)) - exp (LF(p, q)) 



EfS;:::;^ iF ^nU (n?=i^) (li?, % 



(7) 



where ^ denotes that we are summing only over the non-diagonal terms, 
those terms where F(p, q) ^ 0. 

Recall that, without loss of generality, we assumed Ai > ■ • • > Afc. Assume 

qi,i ■ ■ ■ qi, ni 



log 



#0 



Since < exp (35^) and q lf i < exp [Jfl^J, we have 



(8) 



exp N Al 



log 



• • • 9l,m 



> 2 ex P ( ( 1 



40/c 



)iV Al 



• • -Pl,mi. 

which follows from the fact that if m, n are positive integers, and m 7^ n, 
then |log(m/n)| > 1/(2 min(m, n)). Furthermore, for j > 1, for pj G 7^ 
and qj G Qj, then for sufficiently large N, 



exp ( JV Aj 



log 



• • • qj, nj 



• • • Pj,mj 

and so we can conclude that if ( 



< exp ( A" A - ) S A > <- exp ( 2.Y' 1; 
holds, 



|F(p,q)| > exp ((l - -L) - ]Texp ( 2 iV 



> exp 1 



1 

20k 



N' 



for sufficiently large N, since Aj < Ai for all j > 1. 

The contribution of such terms to ([7]) is clearly bounded by 
1 X - 1 



exp ((1 - ^)ATA! 



^ / 7 



< exp 



i"*"l-(i-j5i>"*- 



once more using the fact that 

A: 



y _L#Aj < J_ N ; 
^ A()h ~ 20k 



i=2 



for sufficiently large N. 
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Therefore, as N tends to infinity, we see that terms which satisfy (jHJ) 
vanish. Thus, for a non-zero result in the limit, we must have 

log ( ShLllS^L) = o 

\Pl,l ■ ■■Pl,m 1 J 

That is, we must have (pi,qi) £ T>\, where T>\ is defined in ([5]). 
The terms which might possibly contribute to ([7]) are 

exp (i2F(p, q)) - exp (LF(p, q)) 



E 



K fc!pfeV» ^e^nu [u,;.^, 

q2eQ2,— ,<ik£Qk 

The same argument as above, shows that the terms with 



log ( Q2 > 1 ■ ■ ■ g2 '" 2 ) + 

P2,l • • -^2,7712 



contribute 



E 



( P l,q2)G2?i Ydl^lPl/) (II"=l9l^) 

If (pi,q2) G £>i, then ([6]) shows that 

£ , 1 — « (log aT j 

( P i,q 2 )GOi J{Y\T=iPhe) (n™=i ffM) 



which is negligible compared to exp ( — (1 — igj;)-^ 2 )- Hence, finally, these 
terms do not contribute. 

Repeating the argument for Xj, j = 3, 4, . . . , k we see that any term with 



has a vanishing contribution to the large- iV limit. Therefore, the main term 
must come from those terms for which 

• • • qj,nj 



log _^ ±£^_ = o 

\Pj,l---Pj,m j J 

for all j. Such terms are the diagonal terms, and their contribution has been 
calculated above. This completes the proof of Lemma [6j □ 

Proof of Theorem \4\ Recall that 

logC^ + rue^) 

L x (N,u) = -== . 

Vlog N 



10 

and 
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P(X,n;k,N, u) 



E 



p<exp 



Let 



e(A, n) = e(A, n; k, N, u) 

= L x (N,u) - P(X,n;k,N,u) 

so that, if we write T = exp(iV^), then changing variables to t = Tu, 



e(Xj,mj)\ 2kmj du 



1 1 

(logN) km i T J T 



2T 



logC(| + it)- £ 



= O 



VP 



2krrij 



dt 



(kmj) 4km i e Akm J 
(logA^) fcm i 

by Theorem [5j Since the rrij are fixed, this tends to zero as N — > oo. 
Consider 



(9) 



[]L Aj .(iV, U HL A .(iV, U ) nj ' If /'iA r ///,;'" /'!A r ^;"' 



du (10) 



Writing L Aj . (iV, u) in terms of P(Xj,rrij) and e(Xj,rrij), we see that the term 
inside the modulus signs equals 



e n 



0<ai<mi,...,0<aj ; <m,fe j- 
0</3 1 <m,...,0<P k <n k 



/'(A,., /,,;"'< ">< (Xj. '/'/)" //,) J ((A,.//,) 



The integral of this in ()10p is clearly bounded by 



E 



0<cri <mi ,..-,0<afc<r?i£, 
0</3i<ni,...,0</3 fe <n fe V 

2 k 



n 



mi \ / n 



nj-0i 



]I\P(\ j ,m j )\ m '- a i\e(\ j ,m j )\ a ' P(\j, nj ) ' ' e(X 3 , nj ) 



1 3=1 



du 
(11) 
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A version of the generalized Holder inequality states that 
r2 k 



[ n 1^11^1^,11^1 du 

.7 = 1 



3=1 



\Aj\ 2kr i du 



l/(2kr s ) 



\Bj\ 2ks i du 



l/(2ksj) 



\Cj\ 2kt i du 



l/(2fctj) 



\Dj\ 2ku i du 



l/(2kuj) 



so long as — + — = 1 and ±- + — = 1 for all j = 1, . . . , k. 

Choosing rj = rrij/(mj — ay) and Sj = rrij/ctj, and tj = rij/(nj — f3j) and 
Uj = Uj/Pj, we see that we may bound the above integral by 



n(^V(A„^)i 2 ^ d U 



Tfh~ — Oca 



\e(Xj , mj)\ 2kmj du 







J ' J 

2krij \ 2knj i 






2knj 


t 


P(Xj,rij) 


du) ( 


I 


e(Xj,nj) 





du 



Prom ©, if atj ± 0, 



lim 



and from Lemma [6] we have 

r-2 



leiXj^j)] 2 ^ du 



m.A —a a 



P(Xj , TTlj ) 1 2km i du 



<C 1 



Since the sum in (fTl]) is over those ay and (3j such that ^ otj + (3j > 1, there 
must be at least one j with a non-zero ctj or (3j. Hence, as TV" — > oo, all the 
terms in (jlip tend to zero. The sum is over a finite number of terms, so we 
may conclude that 



lim 



H L\. (N, u)^L Xj (N,uf ] - J] P(\ j ,m j ) m *P(\ j ,n j ) ni 

3=1 3=1 

which implies 

2 k k 



du = 



lim / f[L x .(N,u) m 'L Xj (N,u) ni du = lim [[Pt^m/iPtA^/ 1 

assuming the limits make sense. Therefore Theorem|l]follows from Lemma[6l 

□ 
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